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ABSTRACT 

The  generation  of  gravity  waves  in  a  frontal  zone  is  investigated 
with  a  numerical  model.   The  solutions  show  that  gravity  waves  can  be 
generated  in  a  frontal  zone  by  the  application  of  a  heating  function 
which  varies  sinusoidaly  with  time,  and  that  such  waves  propagate  both 
ahead  of  and  behind  the  frontal  zone.   The  amplitude  of  the  waves  are 
directly  proportional  to  the  magnitude  of  the  forcing  function  and  to 
the  period  of  sinusoidal  oscillation.   Even  with  large  heating  rates,  no 
hydraulic  jumps  are  formed. 
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I.  • INTRODUCTION 

Meteorologists  involved  in  numerical  weather  prediction  have  for  some 
time  been  filtering  gravity  waves  out  of  their  prediction  models.   This 
has  been  done  because  high  frequency  inertial-gravity  oscillations  could 
be  generated  which  would  give  unrealistic  pressure  changes,  obscuring  the 
meteorological  forecasts.   Gravity  waves  may  be  meteorologically  signifi- 
cant however  in  explaining  irregularities  in  the  ionosphere,  and  in  some 
mesoscale  phenomenon  in  the  troposphere.   Hines  (1960)  has  suggested  that 
internal  atmospheric  gravity  waves  may  account  for  many  of  the  irregulari- 
ties observed  in  the  lower  ionosphere  and  that  these  waves  may  have  their 
origin  in  the  large  energy  regions  of  the  lower  troposphere.   Gossard, 
Jensen  and  Richter  (1971)  have  studied  gravity  waves  of  small  amplitude 
with  a  high-resolution  FM/CW  radar  sounder.   Breeding  (1972)  investigated 
the  creation  of  unstable  regions  in  the  atmosphere  by  internal  gravity 
waves.   Gossard  (1962)  has  listed  the  following  three  mechanisms  for 
generating  internal  waves  in  the  troposphere: 

"1.  Internal  waves  can  be  generated  as  standing  waves  in  the  lee  of 
topographic  features. 

"2.  Gravity  waves  are  often  caused  by  a  traveling  boundary  between 
two  types  of  air  of  different  density.   This  type  of  tropospheric  wave 
train  is  short,  lasting  less  than  two  hours.   It  is  a  fair-weather 
phenomenon,  rarely  associated  with  storm  fronts. 

"3.  Oscillations  which  last  for  many  hours  are  generally  associated 
with  tropospheric  storms  and  synoptic  scale  features.   These  oscillations 
of  long  duration  are  usually  generated  by  storms  and  frontal  disturl  nces 
but  the  presence  of  a  stable  layer  in  the  lower  atmosphere  is  requi.vid." 
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Tepper  (1952)  postulated  that  gravity  waves  may  be  generated  in  the 
vicinity  of  a  moving  frontal  zone  if  the  zone  moves  into  an  area  where 
a  stable  layer  exists.   As  the  waves  move  out  they  steepen  into  hydraulic 
jumps.   He  felt  that  large  vertical  motions  caused  by  the  jumps  could 
lead  to  severe  weather  such  as  is  found  near  squall  lines.   Kurth  (1971) 
studied  hydraulic  jumps  in  a  rotating,  continuously  stratified  atmosphere. 
Numerical  simulation  of  squall  line  formation  was  carried  out  by  Sasaki 
(1959)  and  Ogura  and  Charney  (1962).   One  of  the  objectives  of  the  present 
study  was  to  determine  if  such  effects  could  occur  near  a  properly  simulated 
frontal  zone. 

Since  the  advent  of  satellite  pictures  and  the  technique  of  producing 
film  loops  showing  "weather-in-motion"  it  has  become  evident  that  high 
speed  gravity  waves  move  out  from  the  face  of  frontal  zones  (Fujita  and 
Bohan,  1967).  '   * 

In  this  study,  the  non-geostrophic  frontogenesis  model  of  Williams 
(1972a)  is  modified  in  an  attempt  to  generate  gravity  waves  in  a  frontal 
zone. 

In  section  II,  the  Boussinesq  equations  are  given  and  the  modeling 
relations  are  introduced.   In  the  model  the  time-dependent  quantities 
are  functions  of  y  and  z  only.   The  basic  wind  deformation  field  is 
independent  of  time  and  height.   The  numerical  procedures,  boundary 
conditions  and  initial  conditions  are  given  in  section  III.   Numerical 
results  are  presented  and  discussed  in  section  IV.   Conclusions  and 
recommendations  for  further  investigations  are  given  in  section  V. 
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II..  BASIC  EQUATIONS 

In  this  experiment  the  Boussinesq  equations  are  employed  and  the 
domain  is  bounded  by  two  rigid  horizontal  planes.   The  compressibility 
of  the  atmosphere,  which  is  neglected  in  the  Boussinesq  approximation, 
should  not  be  of  qualitative  importance  since  the  density  scale  height 
in  the  atmosphere  is  much  larger  than  the  thickness  of  typical  frontal 
zones.   The  replacement  of  the  tropopause  by  a  rigid  horizontal  surface 
will  induce  large  errors  in  this  region  and  prevents  gravity  waves  from 
propagating  out  of  the  troposphere,  but  the  resulting  errors  near  the 
lower  boundary  should  be  small.   Basic  equation  development  follows  the 
work  of  Williams  (1972a). 

The  hydrostatic  Boussinesq  equations  can  be  written  in  the  following 
form  when  the  earth's  rotation  and  horizontal  eddy  viscosity  are  included, 

4f-+  v-(©V)  t  ^^-/^ze  +q,     (11.2) 

V-V  +  -ff7  =  O  ,     (II. 3) 
"|f    -     3©/@, ,   (H.4) 

where  0  is  a  constant  reference  potential  temperature;  p^  a  constant 
reference  pressure,   k  =  R/(Cp)  ;  @=   T(po/p)   -  0O  ,  the  departure  of 

the  potential  temperature  from  Q0 ,  0=  Cp  @0  (p/po)  +  gz,  the  pressure 

2  .!_ 

function,  A  \7  V,  the  horizontal  diffusion  of  momentum,  A_\7  0,    the 
m  9   v 

horizontal  diffusion  of  heat.,  and  Q  the  heating  function. 
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An  exact  steady  solution  to  this  set  of  equations  is  given  by 


V  =  D  (xi-yj) 

w  =  0 

^  =  J  =  ["_D2(x2+y2)/  2    -f  Dxy 

0-=   0 


,  (II. 5) 


where  D  is  a  constant.   Cartesian  coordinates  have  been  introduced 
and  f  is  taken  to  be  constant.   The  horizontal  velocity  is  a  field  of 
pure  deformation  and  this  deformation  is  given  by  2D. 

If  departures  from  the  fields  (II. 5)  are  independent  of  x,  then  they 

will  remain  independent  of  x.   Thus,  the  dependent  variables  are  subdivided 

as  follows: 

V  =[d  x+  u(y,z,t)J  i+  [-Dy+v(y,z,t)J  j 

w  =  w(y,z,t) 

6>=  e(y,z,t)  | 

p=    $>(x,y)  +  7r(y,  z,  t) 

If  these  expressions  are  substituted  into  Eqs.  (II. 1)-(II.4) ,  they 

become: 

m  +-  lt£  ±.  H*^)    >vu  +  Vi^  -  sir  -  n  A^ 

%+  &  =° 

^^       ©• 

where  V  =  -Dy.   The  boundary  conditions  are 

.  w  =  0,    z  =  0,  H  ,   (11.12) 

where  H  is  the  distance  between  the  horizontal  plates. 

If  the  vertical  average  of  a  quantity  is  defined  as 

W 


,    (II. 

7) 

,      (II. 8) 

,      (II 

9) 

,      (II 

10) 

,       (II. 

11) 

<(  ))  -!/( )^> 
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and  the  hydrostatic  equation  (11.11)  is  integrated  with  respect  to  z 
and  the  vertical  mean  is  removed,  then 

*-<«>=  %[{&* -</&*>]      .(n.13) 

The  y  equation  of  motion  may  now  be  written  in  terms  of   TT  —  \*|T/  • 
If  the  vertical  average  of  (11.10)  and  the  boundary  conditions  (11.12) 
are  used,  then 

by       N 

This  equation  states  that  the  total  disturbance  y  mass  flux  is  independent 
of  y.   If  the  other  variables  have  proper  symmetry  it  follows  that  \xr  / 
must  be  an  odd  function  of  y  which  leads  to 

<V)  =  °  ■  (H.14) 

If  the  vertical  average  of  (II. 8)  is  taken,  and  (11.12)  and  (11.14) 
are  used ,  then 

i<yy>  +  1<2?1  +  f  <;u>  =  o  m 

^  by  (11.15) 

This  result  is  subtracted  from  (II. 8)  which  yields 


*-v- 


Eqs.  (II. 7),  (II. 9),  (11.10),  (11.13)  and  (11.16)  form  a  complete  set 
which  can  be  solved  by  a  pure  marching  process. 


12 


III.   NUMERICAL  PROCEDURE 

The  arrangement  of  variables  and  the  finite  difference  approximations 
are  the  same  as  those  used  by  Williams  (1967).   In  order  to  close  the 
problem  computational  boundaries  must  be  introduced  in  y.   Since  the 
disturbance  velocity  should  die  out  at  a  sufficient  distance  from  the 
axis  of  dilatation,  then 

\r(  +  Y,  z,  t)  -  0  .  (IH.l) 

However,  there  is  appreciable  inflow  across  these  computational  boundaries 
since  V(+  Y)  =  +  DY.   The  quantities  u  and  0  which  are  advected  across 
the  boundaries  must  be  specified  independent  of  the  interior  values  if 
computational  stability  is  to  be  maintained  (Platzman,  1954).   Thus,  the 
following  boundary  conditions  are  used: 

u[+  (Y+^y/2. 

6[±   (Y+Ay/2, 


,z,t]  =  u[+(Y+Ay/2),  z,  o]    ) 

>   .      (III. 2) 
,z,t]  =  0[+(Y+Ay/2),  z,e]    J 


The  computational  boundaries  y  =  +  Y  are  placed  between  grid  points 
so  that  the  above  conditions  are  actually  applied  at  y=+  (Y+Ay/2). 
These  conditions  are  replaced  by  time  specified  functions  in  some 
experiments. 

These  boundary  conditions  developed  a  weak  nonlinear  instability 
near  the  boundaries  but  this  was  removed  by  introducing  a  forward  step 
every  72  time  steps. 

The  initial  temperature  field,  from  Williams  (1972a),  is  given  by 
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were  c*T  -ftf  H~' ($&,,  d6x  /  ^  2r  )  •   The  quantity  <)%  /^ 

which  is  constant,  is  the  initial  stability  and  A  is  one-half  the  total 
horizontal  temperature  variation.   The  initial  x  component  of  the  velocity 
is  given  by 

«(•*,»,*)  »&^C*-HA)«cAC~?).   (in.*, 

The  initial  afield  is  obtained  from  the  quasigeostrophic  equation. 
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IV.   NUMERICAL  SOLUTIONS 

The  finite-difference  forms  of  Eqs.  (II. 7),  (II. 9),  (11.10),  (11.13) 
and  (11.16)  are  solved  by  a  pure  marching  process.   Experiments  are 
conducted  by  varying  boundary  conditions  and  heating  functions.   The 
Williams  model  is  allowed  to  run  for  four  days  by  which  time  a  frontal 
zone  has  stabilized  as  described  by  Williams  (1972b).   At  four  days, 
boundary  conditions  are  changed  or  heating  is  introduced  and  the  model 
runs  for  an  additional  four  days.   The  parameters  which  describe  the 
various  experiments  are  given  in  Table  I. 

The  procedure  for  initializing  the  heating  function  G?  is  as  follows: 

At  the  end  of  four  days,  the  vertical  velocity  field  is  scanned  at 

z=2.5  km  (level  three)  to  determine  the  area  of  maximum  vertical  velocity, 

which  corresponds  to  the  leading  edge  of  the  frontal  zone.   A  region  R 

is  defined  which  is  the  region  where  the  heating  function  is  applied. 

This  region  varies  in  width  and  vertical  extent  depending  on  the 

experiment.   The  function  is  applied  as  follows: 

0   =  C  g(t)  G(y,z)  in  region  R, 
W  =  0  elsewhere. 

Here  C  is  the  maximum  heating  added  and  g(t)  takes  the  value  1.0,  or 

the  following  sinusoidal  variation 

g(t)  =  s  in  (cot) 

where  c^>  is  the  frequency.   The  function  G(y,z)  takes  either  the  value 

1.0  for  block  heating  or  the  value  determined  by  the  sine  squared  function 

r(      .  .   2  (      ir[y-  %M  ] ( 

G  (y ,  z )  =  s  in  1      ■  -  ■  \    ■  v  -.  ■ : y 

J  width  of  heating  zone    f 

In  later . experiments,  a  cooling  function  is  applied  to  the  upper  le  3ls 

of  the  model  to  simulate  evaporation.   This  cooling  is  only  activated 
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TABLE  I.   Variation  of  parameters 
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when  there  is  heating  and  it  is  zero  otherwise.   The  magnitude  of  the 
cooling  is  equal  to  20%  of  the  heating  below. 

Run  number  one  is  a  control  run  with  no  variation  of  parameters. 
Maximum  vertical  velocity  and  the  approximate  position  of  the  frontal 
zone  is  shown  in  Fig.  la.   Departure  potential  temperature,  which  is 
0-  Q-   c^-j-  where  d^x.  /^  2:   is  the  initial  stability,  is  shown  in 
Fig.  lb.   The  divergent  part  of  the  wind  at  four  y  grid  points  is 
shown  in  Fig.  2a  and  Fig.  2b.   Fig.  3a  and  Fig.  3b  show  the  variation 
in  the  departure  potential  temperature  in  the  control  run  and  the  same 
quantity  in  run  number  three  where  a  stable  layer  is  introduced  at  the 
beginning.   In  both  cases,  all  waves  are  virtually  damped  out  at  the  end 
of  eight  days. 

Run  number  two  tests  Tepper's  hypothesis  that  gravity  waves  will  be 
generated  if  a  frontal  zone  moves  into  an  area  where  a  stable  layer  exists. 
A  stable  layer  is  introduced  at  the  end  of  four  days  ahead  of  the  frontal 
zone  at  the  -Y  boundary  between  z=1.5  km  and  z=2.5  km  where  °     1/dZ 
was  set  equal  to  8°K  km   .   No  gravity  waves  or  jump  zones  were  noted. 

Run  number  three  is  the  same  as  the  control  run  except  for  the 
introduction  of  a  stable  layer  at  both  the  +Y  boundaries  at  the  beginning 
of  the  program.   As  in  run  number  two,  the  layer  is  inserted  between 
z=1.5  km  and  z=2.5  km  and  ^i/)i-    is  set  equal  to  8°K  km   between 
these  two  levels.   By  the  end  of  four  days  this  stable  layer  had 
propagated  throughout  the  model  between  levels  z=1.5  km  and  z=2.5  km. 
The  presence  of  a  stable  layer  throughout  the  model  did  not  appear  to 
generate  gravity  waves  or  create  jump  zones. 
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Fig.  la.   Vertical  velocity  w  maxima  and  the  position  of  the 
front  ( ).   Run  number  1. 
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Fig.  lb.   Departure  potential  temperature  69'  at  level  z=1.5  Ian. 
Run  number  1. 
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Fig.  2a.   Divergent  part  of  the  wind  v  as  a  function  of  time  at  four 
levels:  z=.5  km  (1),  z=2.5  km  (3),  z=5.5  km  (6)  and  z=8.5  1cm  (9) 
Run  number  1. 
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Fig.  2b.   Divergent  part  of  the  wind  v  as  a  function  of  time  at  four 
levels:  z=.5  km  (1),  z=2.5  km  (3),  z=5.5  km  (6)  and  z=8.5  km  (9). 
Run  number  1. 
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Fig.  3a.   Departure  potential  temperature  &  '  as  a  function  of  tin;. 
four  levels:  z= .  5  km  (1),  z=2.5  km  (3),  z=5.5  km  (6)  and  z=8.5  km  (9) 
Runs  number  1  and  3. 
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In  run  number  four,  a  heating  function  of  10°K  day   is  introduced  in 
block  form  and  is  120  km  wide,  with  vertical  extent  from  z=.5  km  to 
z=4.5  km,  and  it  does  not  vary  with  time.   This  roughly  approximated  the 
release  of  latent  heat  which  would  be  typical  of  convective  activity  in 
a  frontal  zone.   As  can  be  seen  in  Fig.  4a  and  Fig.  4b  an  increase  in 
vertical  velocity  and  a  definite  wave  pattern  in  the  departure  potential 
temperature  &  '  ahead  of  the  frontal  zone  was  noted.   The  gravity  waves 
appeared  to  be  damped  as  the  distance  from  the  frontal  zone  increased. 
An  additional  result  of  this  experiment  was  that  the  heating  in  the  warm 
air  greatly  strengthened  the  temperature  gradient  in  the  frontal  zone. 
This  coincides  with  the  work  of  Eliassen  (1959)  concerning  the  formation 
of  fronts  in  the  atmosphere. 

In  run  number  five  the  heating  function  took  the  shape  of  a  sine 
squared  curve  with  maximum  heating  occurring  in  the  center  of  a  300  km 
wide  zone  with  vertical  extent  from  z=.5  km  to  z=4.5  km.   This  gave  a 
smoother  heating  function  in  the  frontal  zone.   The  results  were  similar 
to  run  number  four  except  that  the  vertical  velocities  and  gravity  wave 
amplitudes  were  reduced  slightly,  as  might  be  expected. 

In  run  number  six  the  heating  function  was  allowed  to  vary  sinu- 
soidaly  with  time  which  approximates  the  build-up  and  decay  of  convective 
activity  in  the  vicinity  of  a  frontal  zone.   The  frequency  for  a  gravity 
wave  with  a  wave  length  equal  to  the  scale  of  the  heating  zone  gave  a 
period  of  9.9  hours.   As  a  result  of  this  variation,  vertical  velocity 
maxima  appeared  at  Y  2? +1300  km  which  can  be  seen  in  Fig.  5a.   Other 
bands  of  vertical  velocity  relative  maxima  appear  but  only  the  maximum 
values  are  shown  in  this  and  subsequent  figures.   Departure  potential 
temperature  fluctuations  were  also  noted  at  YS/+1300  km  as  can  be    ^n 
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Fig.  4a.   Vertical  velocity  w  maxima  and  the  position  of  the 
front  ( ).  Run  number  4. 
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Fig.  4b.   Departure  potential  temperature  9  '  at  level 
z=1.5  km.   Run  number  4. 
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in  Fig.  5b.  More  significantly,  the  gravity  waves  generated  in  the 
frontal  zone  propagated  in  both  directions  from  the  front  and  extend  to 
the  boundaries  of  the  model.   This  can  be  seen  in  Fib.  6  which  shows 
definite  wave  propagation  in  the  departure  potential  temperature  field 
at  various  locations  on  the  Y-axis. 

In  run  number  seven  the  heating  function  was  reduced  to  3°K  day 
in  order  to  determine  if  the  patterns  found  in  run  number  six  would  still 
be  evident.   Gravity  wave  generation  and  propagation  still  occurred 
but  the  magnitude  of  the  waves  was  greatly  reduced.  Vertical  velocity 
values  also  decreased  and  the  maximum  values  returned  to  the  vicinity 
of  the  frontal  zone.   This  also  was  the  first  run  where  the  heating 
function  was  raised  and  extended  from  z=2.5  km  to  z=6.5  km.   This  slight 
modification  was  made  to  account  for  the  fact  that  convective  activity 
in  cumulus  clouds  does  not  extend  to  the  surface. 

In  run  number  eight,  the  heating  function  was  increased  to  30°K  day  ' 
in  order  to  increase  the  nonlinear  effects  in  the  gravity  waves.   Vertical 
velocity  maxima  remained  at  Y^+1300  km  but  the  magnitude  of  the  vertical 
velocities  were  greatly  increased.   Also  the  first  evidence  of  the 
formation  of  standing  waves  was  observed.   Well  defined  bands  of  vertical 
motion  maxima  were  also  noted. 

The  heating  function  was  increased  to  50°K  day""'  in  run  number  nine 
and  as  can  be  seen  in  Fig.  7a  and  Fig.  7b,  regions  of  strong  vertical 
motion  were  noted  throughout  the  model  as  well  as  corresponding  variations 
in  potential  temperature.   Standing  long  waves  appear  to  form  with  a 
wave  length  of  approximately  360  km. 
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Fig.  5b.   Departure  potential  temperature  O  '  at  level 
z=1.5  km.   Run  number  6. 
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Fig.  6.  Departure  potential  temperature  £ T  as  a  function  of  time  ; 
four  levels:  z= .  5  km  (1),  z=2.5  km  (3),  z=5.5  Ian  (6)  and  z=8.5  Ion  (! 
Runs  number  6,  12  (mod)  and  12. 
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In  order   to   determine  what  effect   a  reduction  of   the  horizontal 
eddy  viscosity  would  have  on  gravity  waves,    in  run  number   ten  Am  and 
A^    were   reduced   to   7,000m2   sec-1   from  20,000  m2   sec-1.      This   reduction 
appeared   to   increase  vertical  velocities    and   led  to   standing  waves   of 
approximately   130  km,    as    can  be  seen  in  Fig.    8a  and  Fig.    8b. 

In  runs   eleven   through   thirteen,    the  heating   function  was   held  at 
30°K  day-1  with  Am  and  A^  held  at   20,000m2   sec-1    ,    and   the  stable  layer 
was    altered.      In   run  number  eleven  the  stable  layer  was    removed  entirely 
and   the   areas   of  maximum  vertical  velocities   shifted   from  Y&+  1300  km 
to  Y£  +  1200  km.      In   run  number   twelve  a  strong  stable   layer  with 

d&x/, c)  *         equal   to   12°K  km-1  was   inserted  between  z=1.5  km  and 
z=2.5  km.      The  maximum  vertical  velocities   moved  almost   to   the   extremities 
of   the  model,    away   from  the   frontal   zone.      Run  number  thirteeen  was    a 
repeat   of   run  number   two  where  Tepper's   assumptions  were   tested,    but  with 
the  inclusion  of    the  heating   function.      Fig.    9a  and  Fig.    9b    show  much 
the  same   results   as   previous   experiments  with   the  exception  that  behind 
the   frontal    zone,   where  no  stable   layer  exists,    the  maximum  w  was   moved 
closer   to   the   frontal   zone.      An   interesting   result  of   run  number   thirteen 
can  be  seen  in  Fig.    10   at  Y=-840km  and  Y=-1140   km  where  a   long  wave 
developed   at   z=8.5  km   (level  nine).      This    agrees  with  Alberty   and 
Van   Sickle    (1969)   who   found  evidence  of   a  mesoscale  wave  between  a 
frontal   zone   and   a  severe  storm. 

Run  number   fourteen  was    another   control   run  which   includes    all   the 
changes   made   to   the  original  model.      It   offers   an  opportunity   to  make 
detailed   comparisons  with   run  number  one.      Table   I   contains   details    of 
parameter   changes   made  between   run   number  one    and   run  number    fourteen. 
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Fig.    8b.      Departure  potential   temperature  &'    at   level   z=1.5  km. 
Run  number   10. 
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Fig.  9b.   Departure  potential  temperature  &%    at  level  z=1.5  km. 
Run  number  13. 
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four  levels  z= . 5  km  (1),  z=2.5  km  (3),  z=5.5  km  (6)  and  z=8.5  km  (S  . 
Run  number  13. 
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Fig.  11a  and  Fig.  lib  show  the  now  common  profile  of  w  and  &  '  distribu- 
tions.  Fig.  12a  and  Fig.  12b  show  the  divergent  part  of  the  wind  v  as 
a  function  of  time  for  four  positions  on  the  Y-axis  ahead  of  the  frontal 
zone.   A  phase  velocity  was  obtained  from  the  propagation  of  the  phase 
lines  of  the  velocity  field.   In  this  case  the  phase  speed  was  approxi- 
mately 25msec-l.   This  can  be  compared  to  Alberty  and  Van  Sickles1 
mesoscale  wave  velocity  of  14m  sec-1  from  observed  data.   Fig.  13  shows 
an  interesting  phenomenon  which  was  present  in  all  the  cases  where 
gravity  waves  propagate  away  from  the  frontal  zone.   As  the  waves  at 
level  three  (2.5  km)  and  level  six  (5.5  km)  propagate,  the  B  '  traces  at 
these  two  levels  diverge.   At  the  point  in  time  where  they  first  begin 
to  diverge,  the  0  '  trace  at  level  nine  (8.5  km)  decreases  while  the  £?  ' 
trace  at  level  one  (0.5  km)  increases  slightly.  When  the  69  '  traces  at 
level  three  and  level  six  begin  to  converge  again,  a  long  wave  appears 
to  be  generated  at  level  nine,  with  a  marked  increase  in  £?' .   The  €3' 
traces  at  Y=-1140  km  seem  to  suggest  that  this  divergent/ convergent 
pattern  repeats  itself.   No  conclusions  will  be  drawn  concerning  this 
phenomenon  at  this  time,  but  it  should  be  an  interesting  area  for  future 
theory  and  study. 

In  runs  number  15  and  16  the  frequency  with  which  the  heating 
function  varies  with  time  is  altered.   In  run  number  15  a  frequency  of 
211/52,000  is  used  corresponding  to  a  maximum  every  14.5  hours.   This 
led  to  an  unrealistic  slope  in  the  frontal  zone  and  an  outward  displace- 
ment of  maximum  w,  as  can  be  seen  in  Fig.  14a.   The  departure  potential 
temepratures  exhibited  wide  fluctuations  as  shown  in  Fig.  14b  and 
Fig.  15.   In  run  number  16  a  frequency  of  2 'fl?  20, 000  is  used  which  •  ed 
to  a  wave  period  of  5.5  hours.   The  results  of  this  experiment  v 
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Fig.  11a.   Vertical  velocity  w  maxima  and  the  position  of  the 
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Fig.  lib.   Departure  potential  temperature^'  at  level  z=1.5  km. 
Run  number  14. 
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Fig.  12a.   Divergent  part  of  the  wind  v  as  a  function  of  time  at  four 
levels:  z=0.5  km  (1),  z=2.5  km  (3),  z=5.5  km  (6)  and  z=8.5  km  (9). 
Run  number  14. 
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Fig.  12b.   Divergent  part  of  the  wind  v  as  a  function  of  time  at  four 
levels:  z=0.5  km  (1),  z=2.5  km  (3),  z=5.5  km  (6),  and  z=8.5  km  (9). 
Run  number  14. 
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(9).   Run  number  14. 
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Fig.  14a.   Vertical  velocity  w  maxima  and  the  position  of  the 
front  ( ).   Run  number  15. 
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Fig.  14b.   Departure  potential  temperature  O   '  at  level  2=1.5  km. 
Run  number  15. 
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Fig.  15.   Departure  potential  temperature  0  '  as  a  function  of  time 
four  levels:  z=0 . 5  km  (1),  z=2.5  km  (3),  z=5.5  km  (6)  and  z=8.5  km 
Run  number  15. 
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quite  similar  to  the  results  of  run  number  seven  in  that  gravity  waves 
were  generated  and  did  propagate  but,  even  with  C=30°K  day"    the  waves 
generated  were  of  very  small  amplitude.   The  maximum  vertical  motion  was 
located  in  the  frontal  zone  and  the  departure  potential  temperature 
profile  was  similar  to  that  of  run  number  one. 

It  was  expected  that  if  the  nonlinear  effects  were  important,  that 
the  scale  of  the  waves  would  become  smaller  as  the  waves  moved  out  from 
the  source  region.   In  order  to  further  examine  these  effects,  the 
boundaries  were  moved  from  1800  km  to  2700  km  in  runs  17-20,   These 
experiments  in  fact  do  show  a  decrease  in  scale  as  the  waves  propagate 
out.   If  this  decrease  comes  from  nonlinear  effects  then  it  should  be 
more  pronounced  in  cases  with  larger  heating.   If  run   19  with  C=3°K 
day~l  is  compared  to  run  20  with  C=100°  day~l,  it  can  be  seen  from 
Table  I  that  the  maximum  w  for  run  20  is  less  than  33  tim^s  the  maximum 
w  for  run  19.   If  a  nonlinear  scale  reduction  occurred  in  run  20  then 
the  maximum  w  in  that  experiment  would  be  larger  than  the  maximum  w  in 
run  19,  times  the  ratio  of  the  heating  rates.   This  seems  to  indicate 
that  the  observed  scale  may  not  be  caused  by  the  nonlinear  effects. 
Rather  the  observed  scale  reduction  probably  comes  from  the  contraction 
in  the  basic  deformation  field. 

A  comparison  was  also  made  with  the  results  obtained  by  Kurth  (1971) 
who  found  that  in  a  rotating,  continuously  stratified  atmosphere, 
hydraulic  jumps  will  form  if  a  certain  relationship  exists  between  the 
Froude  number  and  the  Rossby  number.   The  Froude  and  Rossby  numbers 
were  determined  for  runs  seven,  nine,  fourteen,  and  twenty,  using  the 
Froude  number  F=V/[D(g  ft^d  ft /dz)  1/2  ],  and  the  Rossby  number  R c  =  V/fL, 
with  the  result  that  all  four  cases  fell  on  the  jump  side  of  the 
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critical  curve  F=AR   ,  where  A  is  a  constant  in  the  range  of  6.5-7.5. 
Even  though  jumps  should  have  occurred,  it  is  evident  from  the  previous 
discussion  that  they  did  not  occur.   The  conclusion  must  be  that  gravity 
waves,  by  themselves,  do  not  cause  hydraulic  jumps  under  the  conditions 
of  this  study. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  generation  of  gravity  waves  in  a  frontal  zone  has  been  investi- 
gated numerically.   It  was  determined  that  the  presence  of  a  stable 
layer  ahead  of  a  frontal  zone  was  not  sufficient  to  generate  gravity 
waves  or  hydraulic  jumps.   It  was  found  that  the  addition  of  a  heating 
function  on  the  face  of  the  front  strengthens  the  frontal  zone  and 
generates  gravity  waves,  but  that  these  waves  do  not  propagate  away 
from  the  frontal  zone.   It  was  found  that  if  the  heating  function 
applied  to  the  frontal  face  is  allowed  to  vary  sinusoidaly  with  time, 
that  gravity  waves  will  be  generated  at  the  front  and  will  propagate 
in  both  directions  from  the  frontal  zone  and  induce  vertical  motion 
maxima  both  ahead  of  and  behind  the  front.   It  is  believed  that  such 
motions  could  cause  connective  activity  and  condensation  on  the  less 
stable  warm  air  side.   The  amplitude  of  the  gravity  waves  and  the 
magnitude  of  the  vertical  motion  maxima  varied  directly  with  the  heating 
function  and  the  period  of  oscillation.   A  stable  layer  need  not  be 
present  for  gravity  waves  to  generate  and  propagate,  although  the 
existence  and  strength  of  the  stable  layer  appears  to  affect  the 
distance  the  vertical  velocity  maxima  are  displaced  from  the  front. 
The  vertical  velocity  maxima  fluctuated  with  the  scale  of  the  model 
but  these  fluctuations  appear  to  be  caused  by  the  contraction  in  the 
basic  deformation  field.   The  amplitude  of  the  gravity  waves  generated 
did  not  appear  to  be  large  enough  to  induce  nonlinear  effects  which 
might  lead  to  the  formation  of  hydraulic  jumps  or  squall  lines.   Squall 
lines  do  not  appear  to  be  formed  by  a  simple  Tepper  mechanism  but   ay 
still  be  caused  by  topographic  effects. 
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In  future  studies  using  Williams  model,  a  smaller  grid  scale  should 
be  used  in  order  to  obtain  the  fine  detail  which  may  be  missed  using 
a  course  grid.   The  model  would  be  improved  if  moisture  were  parameter- 
ized or  if  it  were  possible  to  fix  the  forcing  function  based  on 
observational  data.   The  latter  is  not  yet  possible  due  to  inability 
to  measure  the  heating  distribution  .   The  model  should  be  modified 
by  replacing  the  rigid  plate  upper  boundary  with  a  boundary  which  would 
allow  gravity  waves  to  propagate  out  of  the  troposphere.   It  is  also 
suggested  that  cooling  be  included  both  above  and  below  the  heating 
region. 

Possibilities  for  further  study  using  observational  data  would  be 
to  determine  if  a  correlation  exists  between  stable  layer  strength  and 
displacement  of  severe  storms  from  a  frontal  zone,  and  also  whether 
divergence  and  convergence  of  potential  temperatures  curves  between 
two  intermedieate  levels  in  the  troposphere  leads  to  long  wave  formation 
near  the  tropopause. 
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